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ABSTRACT 

Solar sub-surface fluid topology provides an indirect approach to examine the 
internal characteristics of active regions (ARs). Earlier studies have revealed 
the prevalence of strong flows in the interior of ARs having complex magnetic 
fields. Using the Doppler data obtained by the Global Oscillation Network Group 
(GONG) project for a sample of 74 ARs, we have discovered the presence of 
steep gradients in meridional velocity at depths ranging from 1.5 to 5 Mm in 
flare productive ARs. The sample of these ARs is taken from the Carrington ro- 
tations 1980-2052 covering the period August 2001- January 2007. The gradients 
showed an interesting hemispheric trend of negative (positive) signs in the north- 
ern (southern) hemisphere, i.e., directed toward the equator. We have discovered 
three sheared layers in the depth range of 0-10 Mm, providing an evidence of 
complex flow structures in several ARs. An important inference derived from 
our analysis is that the location of the deepest zero vertical vorticity is correlated 
with the remaining life time of ARs. This new finding may be employed as a tool 
for predicting the life expectancy of an AR. 

Subject headings: Sun: activity — Sun: flares — Sun: helioseismology 
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INTRODUCTION 



Solar surface flows have been studied over the last three decades using mostly the 
surface observations such as the photospheric magnetic tracers and Doppler measurements. 



These tracers revealed that t he magnetized regions rotate faster than t 



medium of field-free plasma (IGolub fc Vaiana 



1978 



Komm et al. 



1993 



le surroundin g 



Howard 



19961). 



Dopp 
1979; 



er measurements have shown poleward mer idional flows at the solar surface ( iDuvall 



Labonte fc Howard 



1982 



Ulrich et al. 



19881 ). More recently, it has become possible 



to study sub-surface structures and flows subseq uent to the advent of helioseismology 



( iDeubner &: Gough 



1984 



Gough fc Toomre 



acoustic modes of oscillations f lUlrich 



1970 



19911). which probes th e solar interior using 



Leibacher &: Stein 



197ll ). He lioseismic studies 



have revealed that sunspots are rather shaf 



2000 



Basu et al. 



2004 



Couvidat et al. 



the stage of evolution and age. 



ow, near-surface phenomena ( iKosovichev et al. 



20061). and t 



Svanda et al. 



leir rotation rate with depth depends on 



(120091 ) have reported that a majority of ARs 



displays sudden decrease in the rotation speed, compatible with dynamic disconnection of 
sunspots from their parental magnetic roots. Local helioseis mology has further revealed the 



sunspots to be the locations of large flows near the surface (IHaber et al. 



2004 



2002 



Braun et al. 



Zhao fc Kosovichev 



20041). 



Earlier studies have found that magnetic helicity of A Rs follows a hemisp 



leric rule, 



i.e., positive (negative) in southern (northern) hemi sphere ( Pevtsov et al 



hemispheric 


trenc 


Komm et al. 


2007) 



is also reported for flows in ARs (IKomm et al. 



2005 



20011). Su ch a 



Zaatri et al. 



2006 



20071 ). Zonal flows exhibit larger amplitudes in southern hemisphere especially 



at higher latitudes, coinciding with larger magnetic activit y in that heniisphe re. This may 



be attributed to the inclination of the solar rotation axis fiBeck fc Giles 



advect poleward at nearly the same rate as the quiet regions (QRs) (IHindman et al.l 120091 ) 



2(]()m. Also. ARs 



and show convergent horizontal flows combined with cyclonic vorticity, counter clockwise in 
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the northern hemisphere (iKomm et aLll2004j : IZhao &: Kosovichevll2004l ). 



Activity relatec 


varia 


^ons are also re 


sorted in solar surface i 


ows ( 


Vorontsov et al. 


2002; 


Basu & Antia 


2003; 


Ambastha et al. 


2004; 


Howe et al. 


2005; 


Maurva et al. 


2009 


)• 



is believed that sheared flows in sub-surface layers cause sunspot motions that may lead 
to uns table magnetic topolo gies, causing reconnection of magnetic field lines required for 



flares. 



Basu fc Antial ( 120031 ) have found that the maximum meridional velocity of surface 



flows is smaller when the Sun is more active. The maximum unsigned zonal and meridio nal 



vorticities of ARs are correlated with the total X-ray flare intensity (iMason et al. 



20061). 



Furthermore, stee p meridional velocity gradients are found in flaring ARs at the depth 



range of 4-5 Mm (lAmbastha et al. 



2004 



Maurya et al. 



20091 ) ■ which decreased after flares. 



Although inte r nal flows in ARs anc 



Komm et al 



2009 



Maurya fc Ambastha 



QRs have been studied (IMason et al. 



2006 



2010l ). an understanding of their distinctive 



characteristics in flare productive ARs as compared to that in dormant ARs requires further 
investigation. In this letter, we report on the properties of internal flows in ARs of varying 
levels of flare productivity and magnetic complexity observed during the Solar Cycle 23, 
including large ARs such as NOAA 10030, 10484, 10486, 10070, etc., termed as super-active 
regions (SARs). These SARs are found to possess distinctly different characteristics of 
flows in their interiors as compared to the less productive or dormant ARs. We have 
addressed important issues on the relationship of flow characteristics with the life time of 
ARs, magnetic and flaring activities and hemispheric trends. 
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THE DATA AND ANALYSIS 



For selecting our sample of the ARs, we have used archived information on ARs and 
solar activity provided by the web-pages of Solar MonitoTQ. We first identified the ARs 
producing flares of X-ray class > Ml.O using GOES database! during Carrington rotations 
1980-2052 of Solar Cycle 23. Then, we short-listed the ARs within the heliocentric location 
±40° to avoid projection effects and selected 74 ARs, both flaring and relatively dormant, 
that were well covered by the GONG network. For each flaring AR, we chose a single data 
set corresponding to the day of maximum flaring activity. (A full list of the selected ARs is 
not provided here because of the space limitation.) 



We examined the sub-sur: 



products provided by GONG ( iHill et al. 



ace flows in th e interior of the ARs using the ring data 



20031). It utilizes 16° x 16° x 1664 min data-cubes. 



where the first two quantities correspond to the spatial dimensions in degrees and the 
third is the time duration in minutes. The three-dimensional Fourier transform of the 
data-cube gives the power spectrum which exhibits a trumpet like structure in the [fc^;, fc^, uj] 
space. Slices of the trumpet at given fre quencies uj render concentric rings in [k^, ky] 



plane corresponding to different p-modes ( IHill 



19881 ). A wave propagating through the 



medium with horizontal flow is advected to increased (decreased) frequency depending 
on its propagation along (opposite) the flow. The frequency shift Aw for acoustic waves. 
Aw = U • K, gives the distortion in the shape of the rings. Using a proper fitting technique, 
this distortion can be estimated and the corresponding flow velocity is determined. As the 
modes of different wavelengths are trapped at different depths, the flow patterns derived at 
the surface are the weighted average over depths. This concept is us ed for derivin g the flow 



beneath the surface using regularized least square (RLS) inversion (iGough 



19851) adopted 



^|http : / /www . solarmonitor . org 



http : //www. Imsal . com/SXT/plot_goes .html 
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for GONG data fICorbard et al.ll2003h 



We can derive vertical vorticity, w^, to describe the circulation in ARs, as follows 



(V X u) 



dx dy 

Here, Ux and Uy are the zonal and meridional components of flow, respectively, and Uz is 
the vertical vorticity, which is a physical measure of vertical twist in the flow. In order 
to ascertain its significance as flow characteristic, we carried out an analysis of ARs with 
corresponding QRs located at the same latitude and time, but different longitudes. We 
found that the relative values of Uy and Uz in A Rs were indeed sign i ficant ly larger than 



in the QRs. 



t has been shown o bservationally ( ISvanda et al. 



2008 



theoretically (iJouve &: Brun 



horizontal flows. 



Svanda et al 



20091 ) as well as 



2009) th at Maxwell stresses play a significant role in modifying 



( 120091 ) have also shown that flow parameters obtained for 
ARs are mainly influenced by their magnetic fields as compared to QRs where the flows 
correspond to supergranular structures. Therefore, it is clear that the flows below ARs are 
influenced "magnetically" that are otherwise driven purely by hydrodynamical processes, 
i.e., Reynolds stress and pressure gradient. 

We examined the relationship of internal flow with magnetic activity of ARs using the 
magnetic index, MI, obtained from SOHO/MDI 96 minutes averaged magnetograms. For 
this, we extracted the area of interest from the full disk magnetograms corresponding to 
the 16°xl6° spatial patch used in the ring analysis. Then MI corresponding to the AR is 
calculated from the averaged absolute values of magnetic fields over the patch. We also 
calculated the flare index, FI, for each AR by multiplying the X-ray flux, obtained from 
GOES, with the flare duration and then summing the contributions from all the flares that 
occurred during the 1664 minutes' period of the data-cube. 

We estimated the remaining life times (r) of ARs from the information provided by 
solar monitor for the Earthside and MDI farside imaging as follows. We followed the 
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signature of an AR using magnetograms from the reference, central time tr of the data cube 
until the time when the AR disappeared. In case when the AR rotated past the west-limb 
of the Sun to the farside, we followed up its signature using the acoustic hologram images 
until its disappearance at time te- The remaining life time (r) of an AR is then estimated 
from the difference [te — tr)- 

3. RESULTS AND DISCUSSIONS 

The main results of our analysis for the sample of ARs are presented in Figures [1]- HI 
Meridional velocities, Uy, obtained for two of the ARs, NOAA 10030 (northern hemisphere) 
and 10486 (southern hemisphere), are shown by the solid curves in Figure [H The velocities 
are obtained from the surface to a depth of 14 Mm. The data sets used for these ARs 
correspond to 15 July 2002 and 27 October 2003, respectively. The dash-dotted curves 
show the greatest rate of change, or the gradient of Uy, along the vertical direction z, i.e., 
duy/dz = d'uy. The magnitudes of Uy and d'uy are evidently more significant as compared 
to the errors in the depth range from 1.5 to 10 Mm. 

The profiles of Uy and d'uy for the two ARs, located in opposite hemispheres, exhibit 
opposite trends with depth. The two extrema near P and Q of My, however, occurred at 
around the same depths for both the ARs. Evidently, d'uy changed sign with increasing 
depth from the surface to a depth of 1.5 Mm where d'uy is positive (negative) in the 
northern (southern) hemisphere. Thereafter, it changed sign again in the depth ranges of 
1.5-5 Mm and 5-10 Mm. From these sign reversals of Uy and d'uy profiles, we can infer 
that there exist three sheared layers in the interior of these two ARs. This appears to be a 
general feature of flare productive and complex ARs of our sample, some more of which are 
shown in Figure [2l 
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We may infer that the depths of 2 and 6 Mm near which the extrema of Uy are located 
correspond to the convective scale sizes of granulation and mesogranulation. For reference, 
these depth levels are marked in Figs. [T]-|5]as "G" and "M", respectively. These depths were 



also referred as the locations of ionization zones of H"*" and He"*" ions by 



(1964) and 



November et al 



( 1l98ll ). However, rece ntly 



Cattaneo et 



( I2OO3I ) have shown this argument to be incorrect. iRieutord et al. 



al. 



Simon fc Leighton 



120011) and 



RastI 



(120001) suggested that 



mesogranu 



2001 



Rast 



ations result from the non-linear interaction of granules (also 



20031) 



Rast 



Cattaneo et al. 



(|2003[ ) suggested that granular scales are determined by a close 



association between the upflows, which sustain radiative loss, and the downflow plumes. 



which initiate the upflow motions. A more detaile d description o: 



scales can be found in the recent review article by 



Nordlund et al. 



surface convection and 



fl2009h . 



Convective cells are parcels of fluid moving from the deeper layers to the upper layers 
and back to the bottom. The vertically rising material must change direction along a 
horizontal flow at two positions, i.e., the bottom and the top of convective cells. In a depth 
distribution of the velocity flow, these turning points may be represented by peaks in the 
velocity profile. The extrema seen in the vertical profiles of Uy, such as the one at "G" 
in Figure (H may be due to instabilities that drive the convective cells. The reason for 
convective driving near these depths is not clear because it depends on many factors, e.g., 
physical properties of flux rope and its ambient medium. 

From our sample we found that 44 ARs displayed two extrema in their Uy profiles while 
the rest, i.e., 30 ARs possessed a single extremum in the depth range 0-10 Mm. A statistical 
analysis of these ARs showed that more complex and flare productive ARs possessed two 
extrema, located at depths of 1.92±0.15 and 4.69±0.30 Mm, respectively. On the other 
hand, relatively dormant ARs possessed a single extremum located at a shallower depth 
of 1.66±0.97 Mm, i.e., with large depth uncertainty. There are some further interesting 
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features associated with the ARs having two extrema as compared to the ARs having only 
one extremum: i) larger mean magnetic field, ii) nearly twice the mean GOES fiux, and iii) 
nearly twice the mean life times. 

The maximum d'uy of ARs in the two hemispheres exhibited a general trend, as evident 
from their distribution with depths (Figure Most ARs lying in the northern (southern) 
hemisphere possessed their maximum negative (positive) gradients in the depth range 1.5-6 
Mm. There are 24 (70%) out of the 34 ARs located in the northern hemisphere having 
negative maximum gradients while 29 (74%) out of the 39 ARs located in the southern 
hemisphere having positive maximum gradients. We also noticed that some of those ARs 
which did not follow this hemispheric trend were new emerging ARs. 

It is clear from Figure ISI^a) that the magnitude of maximum d'uy decreased with depth. 
The fits for ARs in northern (southern) hemisphere are drawn in solid (dotted) line. These 
lines are described by d'uy = — 1.15 + 0. 21z for the northern ARs and d'uy = 1.05 — 0.19^; for 
the southern ARs, where, z represents the depth of the maximum d'uy. We do not clearly 
understand the reason as to why the maximum d'uy should have such a linear relationship 
with depth, as only 25% (21%) of the total variance in d'Uy in the northern (southern) 
hemisphere is explained by the linear regression model. However, the hemispheric trend 
is clearly evident from the averages of d'uy obtained for the northern and southern 
hemispheres, i.e., (-5.87 ±0.72) x lO'S"^ and (5.86 ± 0.70) x 10~'^s~\ respectively. 

The relationship of d'uy with magnetic index (MI) of ARs is shown in Figure [3]^b). 
Although the ARs show a rather large scatter, there still appears to be a correlation between 
these two parameters. The fitted lines drawn for the ARs in the northern and southern 
hemispheres are given by d'uy = —0.22 — 0.31 MI and d'uy = 0.01 + 0.48 MI, respectively. 
These fits reveal that the maximum d'uy of ARs increased with MI. The larger slope of the 
fitted line corresponding to the ARs located in the southern hemisphere indicates larger 
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d'uy for these as compared to the ARs of similar MI located in the northern hemisphere. 



We have computed the vertical component of vorticity, w^, for all the ARs of our 
sample, but show the profiles only for two ARs, NOAA 10030 and 10226 ( Figure HI top 
panel). It is evident that uj^ varies with depth both in sign and magnitude as the flow 
parameters shown in Figure [TJ These profiles also show two extrema of opposite signs at 
the depths near 2 and 6 Mm corresponding to the scales "G" and "M". This implies that 
the flows at these depths are twisted in opposite directions. For NOAA 10030, Uz changed 
direction around the depths ~1.8 and ~12 Mm while for NOAA 10226, the corresponding 
depths are ~2 and ~8 Mm. 

Zero vorticity may be a signature of the flux rope being broken at these depths. It 
is possible that the connection between the magnetic structures at the surfac e and its 



unde rlying roots may get broken through a dynamical disasso ciation process ( 



19941 ). and surface reconnection of the two opposite polarities ( ISchrijver fc Title 



an et al. 



19991 ). But 



th ese mechanisms have been overruled by dynamical disconnection mechanism proposed 



by 



Schiissler fc Rempell ( 120051 ) based upon th e buoyant upfiow of 



Schiissler fc Rempell ( l2005l ) found that the 



asma along the field 



lines. From the simulation of thin flux tubes, 
disconnection takes place at a depth between 2 and 6 Mm, which correspond with the 
depths "G" and "M" marked in Figure ll](a). We suggest that once the flux ropes in the 
interior of an AR are disconnected, its remaining life time would depend upon the depth 
where ujz vanishes, i.e., a;^=0. 

We have determined the depth of zero ujz by a minimization procedure. However, it 
was not possible to accurately approach toward the point of zero vorticity for some of the 
ARs of our sample due to the limitation of the available spatial resolution of GONG data. 
For those ARs which did not show zero Uz in their interior, we took its absolute minimum 
to relate with the remaining life time. Further, in Figure HJ^b), we have plotted the depths 
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of zero vorticity, do, with the remaining hfe time, r. A good hnear trend is evident between 
the parameters do and r from the fitted hne, r = 6.43 + 0.34 , drawn through the boxcar 
averaged points. The fit is found to be good as 72.65% of the total variance in r is explained 
by the linear regression model. 

Figure 111(b) shows that ARs having zero vorticity at deeper levels last longer. The 
Pearson correlation coefficient between the depth and life time of ARs is ~85% for ARs in 
the northern and southern hemispheres. However, for the four SARs, i.e., NOAA 10030, 
10044, 10069, 10486, 10488, the remaining life times are found to be >30 days, i.e., much 
larger than the rest of ARs of our sample. These SARs did not follow a linear relationship 
with depth unlike the others. This may be attributed to a continuous process of emergence 
of new flux tubes in these SARs resulting in enhancing their life. Due to their distinct 
abnormal behaviour, we have excluded these SARs from the fittings. 

We have examined the relationship of flare productivity of ARs with the extrema of Uz- 
It is found that the depth of first extremum is mildly correlated with the integrated X-ray 
flux released during the 1664 minutes' period of the data cube obtained for the ARs. The 
Pearson correlation coefficient calculated between these parameters is 20%. On the other 
hand, no correlation is found between the second extremum of Uz and the integrated X-ray 
flux of the corresponding ARs. 

4. CONCLUSIONS 

From our study of the 74 ARs observed during August 2001- January 2007 of Cycle 23, 
we have discovered statistically significant relationship among the sub-surface flow topology, 
energetics and the remaining life time of ARs. We have found the following important 
results: 
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1. Three sheared layers in the depth range 0-10 Mm were found to exist in 44 ARs 
which revealed the complexity of flow structures beneath the surface of more complex 
ARs. The two extrema in Uy and d'uy profiles of these ARs were found to be located 
at the depths of 1.92±0.15 and 4.69±0.30 Mm. 

2. The ARs having two extrema as compared to the ARs having only a single extremum 
in Uy were found to be more active as they possessed as large as twice the mean 
magnetic field (MI), mean GOES X-ray flux and mean life time. 

3. The extrema of meridional velocity gradients were found to follow a hemispheric 
trend, viz., in the northern hemisphere, 24 (70%) out of 34 ARs had negative gradients 
while in the southern hemisphere, 29 (74%) out of 39 ARs had positive gradients. 

4. ARs of larger MI possessed steeper gradient in meridional velocity profiles. 

5. Flaring activity of an AR is found to be associated with depth of the first extremum 
of vertical vorticity. 

6. ARs having zero vertical vorticity at deeper layers are expected to last longer. 

In summary, we have discovered a new hemispheric trend involving the meridional 
velocity gradient of sub-surface flows. More importantly, life time of ARs appears to be 
correlated with the depth of the deepest zero vertical vorticity. This inference may be useful 
in predicting the expected hfe time of ARs. 

This work utilizes data obtained by the Global Oscillation Network Group (GONG) 
program and the Solar Oscillations Investigation/Michelson Doppler Imager (SOI/MDI) on 
the Solar and Hehospheric Observatory (SOHO). Solar Monitor Data was used to find the 
correct location of the ARs. Wc would like to thank the referees for their comments and 
important suggestions that helped to improve the manuscript. 
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Meridionol velocity Uy (m s ') 



Fig. 1. — Meridional velocity profiles (solid curves) in the depth range 0-14 Mm for two 
ARs (a) NOAA 10030 (northern hemisphere) and (b) NOAA 10486 (southern hemisphere). 
The dash-dotted curves correspond to the meridional velocity gradients. Error bars show 
the corresponding uncertainties at la level. Horizontal dashed hues at the depths of 2 and 6 
Mm, marked G and M, correspond to the convective scale sizes of granules and mesogranules, 
respectively. 
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Fig. 2. — Meridional velocity profiles (solid curves) in the depth range 0-14 Mm for nine 
ARs of our sample. The dash-dotted curves correspond to the meridional velocity gradient. 
Error bars show the corresponding uncertainties at Icr level. Horizontal dashed lines at the 
depths of 2 and 6 Mm, marked G and M, correspond to the convective scale sizes of granules 
and mesogranules, respectively. 
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Fig. 3. — Distributions of the maximum meridional velocity gradient d'uy for 74 ARs with 
(a) depth z, and (b) magnetic index MI of the AR. Filled (open) circle represent the northern 
(southern) hemispherical ARs. The solid (dotted) fitted lines drawn through these points 
indicate the general hemispheric trend of the ARs. 
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Fig. 4. — (a) Vertical vorticity profiles, with error bars, corresponding to two ARs, NOAA 
10030 (northern hemisphere) and NOAA 10226 (southern hemisphere). Horizontal dashed 
hnes at the depths of 2 and 6 Mm, marked G and M, correspond to the convective scale 
sizes of granules and mesogranules, respectively, (b) Variation of the remaining life time of 
the sample of 74 ARs with the depth of deepest zero vertical vorticity. 



